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Abstract. The discovery of solar-like oscillations in the giant star ^Hya (G7III) 
was reported by Frandsen et al. (2002). Their frequency analysis was very limited 
due to alias problems in the data set (caused by single-site observations) . The extent 
to which the aliasing affected their analysis was unclear due to the unknown damping 
time of the stellar oscillation modes. In this paper we describe a simulator created 
to generate time series of stochastically excited oscillations, which takes as input an 
arbitrary window function and includes both white and non-white noise. We also 
outline a new method to compare a large number of simulated time series with 
an observed time series to determine the damping time, amplitude, and limited 
information on the degree of the stochastically excited modes. For § Hya we find 
the most likely amplitude to be ~ 2 m/s, in good agreement with theory (Houdek 
and Gough, 2002), and the most likely damping time to be ~ 2 days, which is much 
shorter than the theoretical value of 15-20 days calculated by Houdek and Gough 
(2002). 
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1. Introduction 

The recent detection of solar-hke oscillations in the G7 giant star ^ Hya 
(Stello, 2002; Frandsen et al., 2002) promises new interesting prospects 
for asteroseismology in this part of the Hertzsprung-Russell diagram. 
The amplitude spectrum of ^Hya, based on a one-month time series 
measured in velocity, showed a clear excess of power within a broad 
frequency envelope, similar to that seen in the Sun and other solar- 
type stars (Bedding and Kjeldsen, 2003). The envelope was centered at 
~ 90 /uHz, as expected from scaling the acoustic cut-off frequency of the 
Sun (Brown et al., 1991). The highest peak in the amplitude spectrum 
was ^ 1.9 m/s. 
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Following these observational results, Houdek and Gough (2002) 
calculated the theoretical damping rate, 77, and amplitude for ^Hya. 
The amplitude was 2 m/s and, based on rj, we calculated the damping 
time (or mode lifetime) as l/(2-Kr]) ~ 15-20 days. 

The autocorrelation function of the amplitude spectrum of ^ Hya re- 
vealed a characteristic frequency separation of 6.8 /xHz, in good agree- 
ment with the expected large frequency separation of the radial modes 
(i.e. I = 0) (Stello, 2002; Prandsen et al., 2002). However, using the 
frequencies extracted from the observed amplitude spectrum, Stello 
(2002) showed that a unique solution for the large frequency separation 
could not be found, due to aliasing. It seemed most likely, however, that 
the correct value was indeed in the range 6.8-7.0 /xHz. The observed 
frequencies could be explained by purely radial modes, but the presence 
of non-radial modes could not be excluded, again due to aliasing (Stello, 
2002; Frandsen et al., 2002). 

Based on a small sample of simulated time series, Stello (2002) 
showed that the significance of the observed frequencies, due to the 
effect from aliasing, was strongly dependent on the mode damping time 
adopted for ^ Hya. Hence, in order to quantify the alias problems, the 
damping time for ^Hya has to be known. However, as pointed out by 
Stello (2002), the damping time of ^ Hya could not be measured directly 
from the observed amplitude spectrum (by fitting Lorentzian profiles) 
because the power spectrum was too crowded. Neither was it possible 
to use the CLEANed spectrum to measure the damping time directly, 
since the number of frequencies and their position in the spectrum 
are not known. Therefore, extensive simulations would be needed to 
estimate the damping time. Furthermore, it should be noted that the 
preliminary simulations performed by Stello (2002) indicated that the 
damping time could be significantly shorter than predicted by theory 
(Houdek and Gough, 2002). It is therefore important to establish a 
more precise determination of the damping time from the observations 
than was done by Stello (2002). 

In this paper we describe a time-series simulator, outlining the the- 
oretical background and the technique used to simulate stochastically 
excited oscillations. We use this simulator to determine the damping 
time of ^ Hya using a new method, which is based on comparing the 
overall structure of the observed amplitude spectrum with a large sam- 
ple of simulations. The method also gives the amplitude and some 
limited information about the mode degree of the oscillations. 

The outline of the paper is as follows: Section 2 introduces the 
fundamental ideas used for simulating the stellar signal from stochastic 
pulsations, and describes the parameters that the simulator needs as 
input. Section 3 describes how the simulated time series of ^ Hya were 
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constructed and outlines how we decided the input parameters. Based 
on a few results of the simulated data, we introduce in Sect. 4 the 
method used to determine the damping time of ^Hya and present an 
optimum fit to the observations. Finally, we discuss the method and 
our results in Sect. 5, which also includes the conclusions. 



2. Simulator 

The simulator described in this paper uses the same fundamental ideas 
(also described by Chang and Gough (1998)) as the light curve simu- 
lator developed for the MONS and Eddington missions (Kjeldsen and 
Bedding, 1998; De Ridder, 2002; De Ridder et al., 2003). 

2.1. Stochastic excitation model 

The stellar signal, S, as a function of time, t, is modeled by 

Sit) = y: sAt) , (1) 

where each Sjy(t) is a continuously re-excited damped harmonic oscilla- 
tor that represents a single oscillation mode. 

In general, a damped harmonic oscillator without re-excitation can 
be expressed as 

s,,{t) = A sin(27ri^i + ^)e'*/'^ , (2) 

where A, (f>, and d are the amplitude, frequency, phase, and damping 
time. Rather than assigning a constant amplitude to each oscillator, 
as in Eq. 2, we instead simulate the re-excitation and damping as 
a 'kicking' and damping of the amplitude A. The amplitude of each 
mode is kicked independently at a rate characterized by the small time 
step, Atj^ick, between each kick. The independence of the re-excitation 
is established by having different phases, chosen at random, for each 
mode, so that the time for each kick would not be simultaneous for the 
different modes. After n kicks the amplitude. A, assigned to a mode is 

An = e-'^'^^'^^/'^An-i+en, (3) 

where e"^*'^''^'^/'^ is the damping factor and is the nth re-excitation 
kick, taken at random from a Gaussian distribution with zero mean 
and standard deviation a^- In order to vary both amplitude and phase 
in time, we generate the time series as the sum of sine and cosine terms 
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(using ^sin(a; + 7r/4) = A/y/2{sm{x) + cos(x))), which are simulated 
independently. Hence the expression for Si,{t) used in this simulator is 



liUj. Mence tne expression tor Si,[t) used m tn 
t,{t) = An^i sin(27rz/t + 4>i,^i) + An,2 cos(27rz/t + ( 

tarpprrpssivp nrnrpss sVinwn in F,n S is asvmnt. 



(4) 



The autoregressive process shown in Eq. 3 is asymptotically stable 
up to second order (i.e. it does not die out or 'explode' as n — ^ oo), 
provided |e~^*'''=k/'^| ^ i which is always true for physically meaningful 
values (i.e. Ai^ick > and d > 0) (Priestley, 1981). In the asymptotic 
limit, when the process has relaxed, the mean is {An) = 0, and the 
variance, , can be expressed in terms of as 

'^^^ - 1 _ e-2Atkick/d • 

Using the amplitude-scaled discrete version of Parseval's theorem 
(see Kjeldsen and Prandsen, 1992) we have for each term in Eq. 4 

1 ^ 



n=l 

N 



= 2^ E sin(27rz/t„ + 0^) 



n=l 
1 ^ 



n=l 

2 



= (6) 

for a large number of measurements, N, where P{i^) is the power at 
frequency u in the power spectrum. 

Finally, combining Eq. 6 with Eq. 5 gives the value of (T^ required 
to simulate a continuously re-excited damped harmonic oscillator, as 
shown in Eq. 4, that has an average power P{iy) = A{u)^ in the power 
spectrum. The resulting expression for is 



= A^Ati,i,i,/d, (7) 



where the first-order approximation of has been used (which assumes 
Atkick ^ d). We obtained Eq. 7 by dividing by ^/2 to produce the 
correct standard deviation of the resulting amplitude for the combined 
sine and cosine term (cf. Eq. 4). 

We repeat the autoregressive process (Eq. 3) for a number of steps 
that is significantly longer than the characteristic relaxation time, to 
let it stabilize before the actual simulations of the stellar signal be- 
gin. As a rule of thumb, the minimum number of initial steps of the 
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Figure 1. Simulated noise-free time series of a stochastically excited and continu- 
ously damped oscillation at a single frequency (Eq. 4). The input parameters are: 
d = 1.6 days, v = 115.7 /xHz, A = 0.71 m/s, and Atkick = 2min. Top panel: Time 
series of 3000 data points sampled at 14.4 min. Middle panel: Close-up of the upper 
panel. Bottom panel: Amplitude spectrum of the time series. 



autoregressive process (before starting the actual simulations) should 
correspond to at least twice the damping time. 

The simulated signal for each harmonic oscillator described above, 
Siy(t) (Eq. 4), may be viewed as the sum of two vectors in the com- 
plex plane. The main vector with length (amplitude) A is anchored 
at the origin, cycling around it at the frequency of the oscillation. 
The excitation vector, with variable length, is anchored at the tip of 
the main vector. Its variable direction relative to the main vector can 
be separated into the two orthogonal components: phase (orthogonal 
to the main vector) and amplitude (parallel to the main vector). An 
example of such a stochastically excited oscillation at a single frequency 
is shown in Fig. 1. These simulations looks similar to observations of 
the Sun (Leifsen et al., 2001), which supports the theory of continuous 
excitation by convective elements (Goldreich and Keeley, 1977; Houdek 
et al., 1999). 
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2.2. Amplitude interpolation 

The algorithm described in Sect. 2.1 calculates the amplitude (Eq. 3) at 
regular time steps, but these do not necessarily coincide exactly with 
the times of the observations. Since the amplitude varies slowly, we 
simply use the value closest in time to each actual observation. The 
sinusoidal oscillation itself (Eq. 4), on the other hand, is evaluated at 
the exact times of the observations. The result is a fast and reliable 
method for simulating stochastic excited oscillations with an arbitrary 
observational window function. 

2.3. Noise calculation 

To make realistic simulations of stochastically excited oscillations, it 
is very important to be able to include noise (Kjeldsen, 2003). In our 
simulator both white noise and non- white noise are included. 

The white noise is generated by a Gaussian random-number genera- 
tor. For each data point, the white noise is divided by the weight factor 
associated with that data point before it is added to the oscillation 
signal. 

The non-white noise is created by first calculating the Fourier spec- 
trum of a white noise source generated by a random generator. Then 
we multiply it by a function that describes the desired profile of the 
non-white noise. The result is then converted back to the time domain, 
producing the time series of the non-white noise source, which finally 
is added to the oscillation signal. 

2.4. Input parameters 

To summarize, the following input parameters have to be supplied to 
the simulator: 

1. A set of frequencies. 

2. The amplitude corresponding to each frequency. 

3. The damping time corresponding to each frequency. 

4. The time of each observation (the observational window). 

5. The weight of each observation. 

6. A set of parameters defining the noise levels according to the chosen 
noise function. 

7. The time step between each re-excitation kick Atkick- 
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3. Simulating time series of ^ Hya 

We used the simulator described in Sect. 2 to generate a large number 
of time series similar to that obtained for ^ Hya by Frandsen et al. 
(2002). The input parameters for the simulator were as follows: 

1. Since the observations can be explained as purely radial modes 
(Stello, 2002; Frandsen et al., 2002), this simple case has been 
chosen for the main part of the current investigation. The input fre- 
quencies, 1^1, ... , were the radial modes from a pulsation model 
(Stello, 2002), rescaled so that the mean frequency separation was 
6.8 /iHz, in agreement with the observations (Stello, 2002; Frandsen 
et al., 2002). 

2. The relative amplitude of each mode was defined according to an 
envelope (see Fig. 2, solid vertical lines). The envelope was obtained 
by smoothing the observed amplitude spectrum, subtracting the 
noise background and normalizing the peak to unity. In order to 
match simulations with observations the width of the envelope was 
made adjustable. We chose to do this by raising the curve to the xth 
power, where x was a free parameter, which was a convenient way 
of changing the width of the envelope without changing the height. 
The envelope was then scaled vertically by a factor Ampg^g^^g, which 
was a free parameter that represents the oscillation amplitude of 
^ Hya. Figure. 2 (dashed curve) shows the normalized envelope for 
x = 2 and Ampg^^ig = 1 m/s. 




100 

Frequency/ /u.Hz 



200 



Figure 2. Three amplitude envelopes: dashed curve: x — 2 (FWHMin = 48.0 //Hz); 
dashed-dot-dot-dot curve: x = 1.2 (FWHMin = 64.0 /xHz); dotted curve: x = 0.7 
(FWHMin = 81.6 /uHz). The vertical solid lines indicate the mode frequencies. 
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Figure 3. Spectral window of the ^ Hya time series. 

3. The damping time, d, was a free parameter and was the same for 
all modes. This is probably a reasonable approximation, since the 
theoretical damping rate shows a flat plateau covering a fairly broad 
range of frequencies (~ 70-130 /xHz) (Houdek and Gough, 2002) . 

4. The observational window function was exactly the same as for the 
actual observations (Stello, 2002; Prandsen et al., 2002). The ampli- 
tude spectrum of the observational window (the spectral window) 
is shown in Fig. 3. 

5. The weight of each observation was l/cr|, where CTj was the noise 
associated with each observed data point. 

6. The noise was the sum of a white and a non-white noise compo- 
nent, where the latter was described by a linear model. The noise 
was therefore specified by three parameters: the white noise level 
(Noise^hitc)) the slope of the non- white noise (Noisegiope)) and a 
scaling of the non- white noise (Noisegcaie)- 

7. The time step A^kick was set to 2 min to meet the requirement that 
Atkick ^ d for all damping times tested in this investigation. 

Before starting the actual simulations, we repeated the autoregres- 
sive process (Eq. 3) for 120000 steps, corresponding to 167 days, which 
is significantly longer than the characteristic relaxation time for all the 
damping times tested in this investigation, to ensure it had stabilized. 
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4. Method and results 

Examples of simulated amplitude speetra with different damping times 
are shown in Fig. 4. It is clear that amplitude spectra based on a 
short damping time have much more densely packed peaks (top panel) 
than do those with long damping times (bottom panel). This arises 




100 

Frequency//iHz 



Damping time: 12 days 




100 

Frequency/ ^Hz 



100 

Frequency / fj,Uz 




Figure 4- Amplitude spectra based on simulated time series. The input damping 
times for the time scries arc shown in each panel. For infinite damping time (bottom 
panel) the input amplitudes were randomized. 
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Figure 5. Amplitude spectrum (observed) of g Hya. Different parameters character- 
izing the spectrum are indicated (see text page 10). 



because the continuous re-excitation of modes introduces slight shifts 
in the phase as a function of time, which shows up in the amphtude 
spectrum as extra peaks sHghtly offset in frequency. The fact that am- 
phtude spectra based on different damping times display such different 
characteristics, as seen in Fig. 4, suggests the possibility of measuring 
the mode damping time from the overall structure of the amplitude 
spectrum. 

We determined how well the simulated time series reproduced the 
observations by using eight measurable parameters that characterized 
different features of the amplitude spectra. In Fig. 5 we show the 
observed amplitude spectrum of ^ Hya, together with the measurable 
parameters that specify different characteristics of the amplitude spec- 
trum. 

The highest peak (MaXpcak) is the highest amplitude found in the 
spectrum, while MaXsmooth is the height of the smoothed amplitude 
spectrum, where smoothing was done twice with a boxcar filter. The 
widths of the boxcars were 20 /zHz followed by 5 /xHz. The two noise 
levels, Noiseiow and Noischigh are measured as the mean amplitude 
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in the frequency ranges 5-45 /xHz and 145-185 /xHz, respectively. The 
number of detected peaks with amphtudes above a threshold of 1.0 m/s 
(S/N> 3.5), denoted Npcaks) is found by CLEANing the amphtude 
spectrum until the amplitude threshold has been reached. Our CLEAN 
process subtracts one frequency at a time (the one with the highest 
amplitude), but recalculates the amplitude, phase and frequencies of 
the previously subtracted peaks while fixing the frequency of the latest 
extracted peak. In this way, the fit of sinusoids to the time series is 
done simultaneously for all peaks. The mean amplitude is measured 
both over the entire frequency range 0-190 /iHz (Amptot) and in the 
central part that is dominated by the stellar excess of power, 40- 
140/LtHz (Ampcen)- Finally, the width of the excess power (Wenv) is 
measured as the FWHM of the smoothed spectrum after the noise 
has been subtracted. Due to the window function, stellar excess power 
'leaks' into the frequency regions where we measure the noise. The 
subtraction of noise in the determination of the envelope width there- 
fore includes some stellar power, making Wenv smaller than the input 
widths (Fig. 2), which therefore should not be compared. 

We now describe how the eight parameters measured from the simu- 
lated amplitude spectra (Maxpeak, MaXsmooth, Noiseiow, Noisehigh, Npeaks 
Ampccrn Amptot, Wenv) are affected by changes in the input parame- 
ters, namely damping time (d), noise (Noise^hite) Noisegiope) NoisCgcaie)) 
amplitude (Ampg^^jg), envelope width (determined by the exponent x), 
and input frequencies (z/i, . . . , i/„). 

4.1. Damping time 

Figure 6 shows the dependence of each of the eight measured pa- 
rameters on the damping time. Values adopted for the other input 
parameters are given in the figure caption. Each plotted point is the 
mean of 100 simulations with the same input parameters but differ- 
ent random number seeds, and the vertical bars show the rms scatter 
over these 100 simulations. This rms is the intrinsic scatter due to the 
stochastic nature of the oscillations, and is therefore the quantity we 
should use to decide whether the simulations match the observations. 
The parameters measured from the observed amplitude spectrum are 
indicated by horizontal dashed lines. 

All the parameters increase as the damping time gets shorter except 
for MaXpcak) which falls off. The fall in MaXpcak is simply because the 
re-excitation spreads the power over more peaks, giving less power in 
each. The relative change in most of the parameters over the range of 
damping times plotted in Fig. 6 is fairly small (< 10%) and is generally 
less than the rms scatter. These parameters are therefore not very sen- 
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Figure 6. The eight measured parameters as a function of the damping time. The 
data points are based on 7 distinct sets of simulated spectra only differing in their 
damping time. Each set comprised 100 independent simulated time series. Each 
plotted point is the mean of the values found from the 100 simulations within each 
set, and the vertical bars indicate the rms scatter that results from the stochastic 
nature of the simulation. The points are connected by solid lines to guide the eye. The 
dashed lines denote the values measured from the observations. The common input 
parameters for the 7 sets are: Amp^^^i^ = 2.1 m/s, x = 2.0 (FWHMin = 48.0 /aHz), 
Noiscwhite = 0.20 m/s, Noisesiope = 1.5, Noisescaie = 3.0, and 18 input frequencies 
with a mean frequency separation of 6.8 /xHz. 
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sitive measures of the damping time but should still be matched with 
the observed values to constrain the other input parameters. However, 
MaXpcak Ei'iid Npcaks change by roughly 50% in the same range and 
in opposite directions, making them the obvious parameters of choice 
for constraining the mode damping time. The correlation coefficient 
between MaXpeak and Npeaks is only p ~ 0.10-0.15, based on a few 
sets of 100 simulated amplitude spectra, where each set had different 
input parameters. These two parameters can therefore be regarded as 
uncorrelated. 

The 2nd, 3rd, and 4th panels in Fig. 6 all indicate that, for all 
damping times, the amount of power in these simulations is too low. 
The bottom panel shows that this is partly because the width is too 
narrow. The input power can be adjusted by changing the input noise, 
amplitude, envelope width, and the number of frequencies (their mean 
separation). We address each of these possibilities in turn in the next 
sections. 

4.2. Noise 

We investigated how the measured parameters changed as a function 
of input noise. We only show the results of changes in NoisCwhite be- 
cause this parameter gave us all the control we needed to adjust the 
simulations to match the noise level in the observations. The two other 
noise parameters were fixed at the values used for Fig. 6. 

Figure 7 shows that the parameters most sensitive to changes in the 
input noise are Noiseiow and NoisChigh, while the other parameters show 
much smaller changes. The input noise can therefore be regarded as the 
final fine tuning parameter, and has less importance for constraining 
the other input parameters. The panel showing the measured noise will 
therefore be omitted in the following plots. 

4.3. Amplitude 

In Fig. 8 we show the dependence of the measured parameters on the 
input amplitude. It can be seen that all parameters show an increase 
with input amplitudes, which is expected due to the increase in power, 
although the change is very small for Wcnv 

It is clear from Fig. 8 that the input amplitude is constrained by 
MaXpeak) which should not be too high, and by Npeaks ! which should not 
be too low. The difficulty in satisfying both constraints simultaneously 
becomes greater for larger values of the damping time. 
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Figure 7. The effect of changing the white noise. These panels are similar to those 
plotted in Fig. 6, but for three different values of the white noise. Each measured 
parameter is therefore shown as three points for every damping time. The values for 
input noise (from top to bottom in each panel) are: Noisewhito ~ 0.2 m/s (as in Fig. 
6), Noiscwhite = 0.1 m/s, and Noisewwte = 0.0 m/s. All parameters show an decrease 
as Noisewhite decreases. Note that the y-axis has been shifted, but not scaled, relative 
to Fig. 6. For clarity, the vertical bars indicating rms scatter have been omitted. 
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Figure 8. The effect of changing the input amphtude. These panels are similar to 
those plotted in Fig. 6, but for throe different values of input amplitude. Each 
measured parameter is therefore shown as three points for every damping time. The 
input amplitudes (from bottom to top in each panel) are: Amp^^^j^ = 2.1 m/s (as 
in Fig. 6), Ampgj,jjij, = 2.5 m/s, and Amp^^^j^ = 2.8 m/s. All parameters show an 
increase as Ampg;.g^jg increases. 
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Figure 9. The effect of changing the envelope width. These panels are similar to 
those plotted in Fig. 6, but for three different values of envelope width. Each mea- 
sured parameter is therefore shown as three points for every damping time. The 
input widths (from bottom to top in each panel) are: x = 2 (FWHMin = 48.0 ^iHz) 
(as in Fig. 6), x = 1.2 (FWHMin = 64.0 /iHz), and x = 0.7 (FWHMin = 81.6 /uHz). 
All parameters show an increase as the envelope width increases. 



stello_complete.tex; 2/02/2008; 8:44; p. 16 




90 p— 
BO E- 
^ 70 E- 
* 60 §- 
50 i- 



— * 



5 10 15 ao 

Damping time/days 

Figure 10. The effect of changing the number of input frequencies. These panels 
are similar to those plotted in Fig. 6, but for three different values of the frequency 
separation. Each measured parameter is therefore shown as three points for every 
damping time. The frequency separations (from bottom to top in each panel) are: 
6.8 /iHz (as in Fig. 6), 5.0 /xHz, and 3.4 /mHz. All parameters show an increase as the 
frequency separation decreases. The 3.4 /xHz simulations corresponds to having one 
non-radial mode per radial mode positioned halfway between each radial mode. 
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4.4. Envelope width 

Another way of injecting more power into the simulations, enabhng 
more peaks to be detected without having Maxpeak increase signifi- 
cantly, is by making the envelope wider (of. Fig. 2). Figure 9 shows 
the dependence of the measured parameters on the envelope width. 
As expected, an increase in envelope width gives rise to an increase in 
all parameters, due to the power increase. MaXpeak is, however, nearly 
unaffected. Thus, by increasing the width we construct simulations from 
which more peaks are detected (Npeaks increases), without affecting 
MaXpeak) hence making the simulations better fit the observations, 
which is also supported by the bottom panel. Since the input width 
is the only input parameter to induce significant changes in Wcnvj we 
see from the bottom panel that the input width can be fixed at a value 
around SO/nHz. 

4.5. Input frequencies 

Finally, wc tested the dependence of the measured parameters on the 
number of input frequencies. The above examples have included only 
the radial modes, so we have simulated the presence of non-radial 
modes by decreasing the mean frequency separation. Since we used 
the radial modes from a pulsation model (Stello, 2002), they showed 
a small scatter of 0.2 //Hz from perfectly uniform spacing. This scatter 
was kept constant as we reduced the frequency separation and included 
more frequencies. 

The results of the change in frequency separation are shown in Fig. 
10. As expected we see an increase in all measured parameters when 
increasing the number of input frequencies (decreasing the frequency 
separation), due to the increase in the power, although Wenv is nearly 
unaffected. 

4.6. Finding the optimum input parameters 

We find the most likely damping time by examining Maxpeak/Maxsmooth) 
which is relatively insensitive to Ampsj,j^ig, envelope width, input fre- 
quencies, and noise, but not to the damping time (see Fig. 11). Thus 
the properties of the amplitude spectra we measure with the ratio 
Maxpeak/Maxgmooth are mainly determined by the damping time. This 
ratio therefore provides a robust measure of the damping time, although 
it is not very precise (due to scatter). From Fig. 11 we see that the 
most likely damping time, given the observations, is d 2± 2 days. In 
Fig. 12 we show the results of simulations with a set of input param- 
eters optimized for d ^ 2-3 days while having the envelope fixed at 



stello_complete.tex; 2/02/2008; 8:44; p. 18 



Simulating stochastic excited oscillations 19 



Amplitude change 





5.0 










4.5 




















i 


4.0 


























iO 

\ 


3.5 
































3.0 
















a 










































2.5 




















2.0 













Envelope width change 




Frequency separation change 




2.0 - 



White noise change 




Damping lime/days 

Figure 11. The ratio Maxpcaii/Maxsmooth as a function of damping time. The four 
panels show the sensitivity of MciXpeak/MajCsmooth to changes to the different input 
parameters. The values of the input parameters are the same as shown in the former 
figures in Sect. 4. 
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Figure 12. The deviation between observations and simulations relative to the rms 
scatter of the simulations (sigma deviation) of all eight measured parameters as a 
function of input damping time. The input parameters (optimized for d ~ 2-3 days) 
are: Amp^^^i^ = 2.0 m/s, x = 0.7 (FWHMin = 81.6 /xHz), NoisCwhite = 0.10 m/s, 
Noisesiopc = 1.5, Noisesoaie = 3.0, and 18 input frequencies with a mean frequency 
separation of 6.8 pHz. 



FWHMin = 81.6 /iHz (cf. Sect. 4.4). All eight parameters are plotted 
in a single plot by showing the difference between observations and 
simulations relative to the rms scatter in the simulations. We find the 
most likely amplitude of ^Hya to be Ampg^^^jg ~ 2.0 m/s. 

When we optimized the input parameters for a damping time of 
~ 15 days, we got roughly the same result as shown in Fig. 12. Changing 
the amplitude or the number of frequencies does not change Maxpeak 
and Npeaks without also affecting the other parameters and worsening 
the fit (cf. Figs. 8 and 10). The problem with a damping time of ~ 
15 days or longer is that we cannot produce enough high peaks while 
keeping the other measures of power down at the observed level. 

It is possible to include a few extra modes (in addition to the 18 
modes with mean frequency separation 6.8 fiKz) in the simulations and 
still have an acceptable fit with the observations because the change in 
the measured parameters is small, relative to the scatter, provided we 
also reduce the input amplitude slightly (see Figs. 8 and 10). 
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5. Discussion and conclusion 

Wc have, in the previous sections, described a method for generating 
reahstic simulations of stochastically excited oscillations, including an 
arbitrary noise function and an arbitrary window function. Using time 
series produced by our simulator we showed how the damping time, 
amplitude, and other mode properties of the oscillations could be de- 
termined by comparing the overall structure of observed and simulated 
amplitude spectra. The method was applied to the single-site time 
series of radial velocity measurement of the red giant star ^ Hya (Stello, 
2002; Frandsen et al., 2002). 

Due to the stochastic nature of the simulated solar-like oscillations, 
we see large variations in the amplitude spectra when the length of 
the time series is not significantly (> 10 times) longer than the damp- 
ing time. Hence, a large scatter is induced for some of the measured 
parameters that describe the characteristics of the amplitude spectra 
(see Fig. 6) . In the case of Hya, we are therefore not able to exclude 
the damping time of ^ 15-20 days calculated by Houdek and Gough 
(2002). However, based on this single dataset, a shorter damping time 
of only a few days seems much more likely (see Fig. 12 diamonds and 
circles). A clear rejection of a damping time of 15 days or longer would 
require a reduction of the scatter seen in, e.g.. Fig. 11 by at least a 
factor of two, and hence, using single-site observations, a time series of 
~ 150 days (assuming that the measured parameters from the observa- 
tions are unchanged). The optimum fit to the observations, assuming 
purely radial modes, gave a maximum amplitude ~ 2.0 m/s, in good 
agreement with the calculations by Houdek and Gough (2002), and 
damping time ~ 2-3 days. Also, due to scatter, it was not possible to 
exclude the presence of a few extra modes beside the 18 modes with 
mean frequency separation 6.8 //Hz. Since the theoretical calculations 
of the mode inertia (Teixcira, 2002) show that the radial modes should 
be excited to larger amplitudes than the higher degree modes, we con- 
clude that our simulations strongly suggest that the 4 Hya amplitude 
spectrum is dominated by radial modes, but with a possible presence 
of a few higher order modes. 

Below we itemize the limitations of the current investigation and 
discuss their consequences: 

Our approximation of using a damping time that does not vary 
with frequency underestimates the damping time for low frequency 
modes (by creating more peaks) while overestimating the damping 
time at the high frequency end (producing fewer peaks) relative to 
the theoretical damping rate (Houdek and Gough (2002); Fig.l). 
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Hence, due to these opposing effects, this should have httle effect 
on our results. 

We included non-radial modes by assuming that they were excited 
to the same amplitudes as the radial modes in the same frequency 
range. A comprehensive treatment of the non-radial modes should 
include both the relative mode inertia, based on a pulsation model, 
and the spatial response of the observations due to projection 
effects for modes of different degree. The technique presented in 
this paper to measure mode life time is probably not sensitive 
enough to make worthwhile a more comprehensive treatment of 
the non-radial modes. 

We tested for any significant effect from the adopted deviation of 
the frequencies from equal spacing by testing both the scatter of 
0.2 iiHz seen in the pulsation model (as in the examples shown in 
Sect. 4) and the larger scatter (0.6 /xHz) seen when the observed 
frequencies were ordered to match the radial modes (see Stello, 
2002). In these two cases, the eight measured parameters were 
found to be nearly unchanged, with a deviation of less than O.la, 
which is the precision with which we know the mean values based 
on 100 independent simulations. 

Furthermore, we tested a very irregular frequency distribution by 
using the 13 observed peaks with S/N> 3.5 (sec Stello, 2002) 
as input frequencies, neglecting possible contamination from alias 
peaks and noise. The results are similar to those shown in Sect. 4. 

We conclude that the regularity in the input frequencies is not 
important for obtaining our current result. 

The observed velocities of ^ Hya were reduced using one reference 
point per night, which produced a high-pass filtering of the time se- 
ries (Frandsen et al., 2002). It would presumably require somewhat 
different values of the input parameters to match the unfiltered 
amplitude spectrum. We expect this would mostly affect the input 
noise and amplitude, by underestimating them, and to a lesser 
extent the damping time. 

We chose to use the smoothed observed amplitude spectrum as the 
frequency envelope because it represents the actual observations. 
We expect that a Gaussian profile could be used for simplicity with 
little effect on the results. 

We used two successive relatively narrow boxcars to smooth the 
amplitude spectrum in order to obtain Maxsmooth- This method 
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preserved the large-scale structure while removing variations on 

small scales. Increasing the boxcar widths by a factor of 2 or 3 
did not produce any significant change in the difference between 
Maxsmooth of simulations and observations. 

- To test the robustness of our method we applied it to measure- 
ments in velocity of aCen A (Butler et al., 2004). Based on a 
plot similar to Fig. 11, we obtain a damping time of ^ 0-5 days. 
A complete analysis, as shown for ^Hya, has not been applied to 
a Cen A. Since a Cen A is very similar to the Sun, one would expect 
the damping time to be only a few days. Our result is consistent 
with both the solar value (3-4 days) and the value of 1—2 days for 
a Cen A found by Bedding et al. (2004). 

In general, we see from our simulations that when comparing the 
estimated amplitude (either from scaling or theoretical calculations) 
with observations, the amplitude associated with a star from a single 
dataset can vary significantly, especially for stars with mode life times 
that are long compared to the length of the time series (cf. Fig. 6). 

In a future paper, a more detailed frequency analysis of ^ Hya will 
be presented, including an analysis of the scatter of frequencies about a 
uniform distribution, which would test the damping time results given 
in this paper. 
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